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Abstract. We study the effect of electron-electron interaction on a two dimensional (2D) 
disordered lattice. For the case of two electrons the analytical estimates are presented showing 
a transition from localized to delocalized states in a way similar to the Anderson transition in 
3D. The localized phase corresponds to large values of the parameter r s , which is determined 
by the ratio of the Coulomb and Fermi energies. The numerical investigations of the spectral 
statistics P(s) in a system with up to N p = 30 spin polarized electrons show the transition 
from the Poisson to the Wigner-Dyson distribution at a total electron energy independent of 
N p . The relation to experiments on the metal-insulator transition in 2D is also discussed. 
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1 Introduction 

Since 1979 it became clear that non-interacting electrons are always localized in a 2D 
disordered potential jlj . At present this result is firmly confirmed by different analyt- 
ical and numerical studies (see the reviews 0, ||) according to which all one-particle 
eigenstates are exponentially localized in the case of orthogonal or unitary symme- 
try. Delocalization is possible only for the symplectic symmetry which is however not 
realized in the absence of spin-orbit interaction fi EJ. Therefore, the discovery of 
metallic behavior in 2D high mobility samples by Kravchenko et. al. Q] attracted 
a great interest of the solid-state community and pushed forward the question about 
a role of electron-electron interactions in the localized phase. Indeed, in many ex- 
periments of different groups §, |, |, % % g, fjl], Q the ratio of the Coulomb 
interaction energy E ee to the Fermi energy ep is characterized by a large dimension- 
less parameter r s — 1/ y/Tm s a* B ~ E ee /tp, where n s is the electron density in 2D, and 
a* B = h 2 ea/m*e 2 , m*, eo are the effective Bohr radius, electron mass and dielectric 
constant, respectively. Typically in the above experiments the parameter r s is ranged 
in the interval 5-40 showing that these experiments are well outside of the perturba- 
tive regime corresponding to r s < 1. Generally, in these experiments the insulating 
increase of resistivity at low temperature T is observed at large r s values (low charge 
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density n s ), while its metallic decrease at low T is seen at lower r s (higher n s ). In ad- 
dition, experimentally it was found that a sufficiently strong magnetic field parallel to 
the 2D plane destroys the metallic behavior Jl3|, |l4| . These results show an important 
role of the spin degrees of freedom. However, we should note that the recent experi- 
ments Jl5[ with n-SiGe samples demonstrate that the metallic behavior persists even 
in very strong in plane magnetic field when all spins are polarized. Another important 
experimental indication was found for a localized phase with a variable range hopping 
(VRH) Jl(| |T^] . Indeed, the experiments show that a prefactor in the exponential VRH 
resistivity dependence on temperature is phonon independent |l7j . For the theoretical 
explication of this result other physical mechanisms of VRH should be found |ll| and 
it is possible that their origin is related to the electron-electron interaction as it had 
beed discussed long ago in fL9f . 

From the theoretical view point the problem of interaction in the localized phase 
is rather nontrivial. Indeed, due to localization of noninteracting states, the two 
interacting particles (TIP), with a short range repulsive/attractive interaction and 
located on a distance of one-particle localization length l\ from each other, always 
return to their center of mass that enhances enormously their interaction. As a result 
the two particles can propagate coherently on a distance much larger l\. This effect 
has been discussed recently by different groups for a short range interaction |^0|, |2l|, |22] , 
p3| . [2-| [25[| (see also the early paper by Dorokhov for a strongly attractive case |26f). 
However, in the experiments discussed above the charge density is very low and the 
distance between electrons (holes) is much larger than a* B (r s 3> 1). In this situation 
the localized electrons interact via the long range Coulomb interaction. Recently, it has 
been shown |27], ^| that in this situation two interacting electrons can be delocalized 
even when they are separated by a distance R^> l\. Moreover, this delocalization goes 
in a way very similar to the Anderson transition in 3D. In this paper we discuss the 
physical origin of this delocalization (Section II). After that (Section III) we present 
the results of extensive numerical studies of up to N p — 30 electrons (with polarized 
spins) on 2D disordered lattice. In agreement with the first results presented in |2^] we 
establish that the ground state remains localized (nonergodic) but at moderate values 
of parameter r s the electrons become ergodic at very low total energy where their 
spectral statistics demonstrates a transition from the Poisson to the Wigner-Dyson 
distribution. The conclusion is presented in Section IV. 



2 Two interacting electrons in the 2D Anderson model 

Let us consider a system of electrons with polarized spins on a 2D disordered lattice 
(Anderson model). Such polarization can be reached by a strong in plane magnetic 
field but also one can speak about spinless fermions. The system is described by the 
Hamiltonian 

<ij> i i>j % ^ 

Here a\(ai) is the fermion creation (annihilation) operator at site i, the hopping be- 
tween the nearest neighbors is V, the diagonal energies Wi are randomly distributed 
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within the interval [—W/2,W/2] and U is the strength of the Coulomb interaction 
with rij being the distance between electrons at sites In this notation m — aja, 
is the occupation number at site i. The electrons (particles) are moving in a 2D cell 
of size L x L with periodic boundary conditions. The Coulomb interaction is taken 
between electrons in one cell of size L and with 8 charge images in nearby 8 cells as 



in 27 . The number of particles N p and the cell size were varied within the intervals 
2 < N p < 30 and 8 < L < 31. With the notations of Eq. (1), the parameter r s is 
given by r s — U/(2V^/ttv), where v — N p /L 2 is the filling factor and ep — AttvV. 
The majority of our data have been obtained for U/V = 2, 0.0048 < v < 0.03 that 
corresponds to 3.22 < r s < 8.14. When changing N p the filling factor was kept 
approximately constant (nearest rational value) by the appropriate choice of L. 

For U = all states in the Hamiltonian (1) are localized and the one-particle lo- 
calization length varies exponentially with W in the limit of weak disorder: lnii ~ 
(V/W) 2 . The numerical diagonalization j2?], £28| allows to determine the inverse par- 
ticipation ratio (IPR) for one particle £i, which approximately gives the number of 
lattice sites contributing in one eigenstate. The values found at L = 24 for different 
disorder strength at the ground state and the middle of the band are correspondingly: 
£i = 3.4 and 4.2 (W/V = 15); 5.2 and 11.2 (W/V = 10); 8.2 and 36.7 (W/V = 7); 
13.5 and 84.2 (W/V = 5). This shows that the one-particle states at low energy and 
even at the band center (except maybe W/V — 5) are well localized (£i <C L 2 ). 

To understand how two electrons interact in the 2D localized phase we should esti- 
mate the interaction induced transition rate T2 between unperturbed (U = 0) localized 
eigenstates. Following |27| ^8| let us assume that the distance between electrons R is 
much larger than li (R 3> lx) and their energy is in the middle of the band. In this 
case the two-body interaction has a dipole-dipole form since the lower order terms 
only slightly modify the one-particle effective potential. The two-body matrix element 
between noninteracting eigenstates is then U s ~ Ulf/R 3 x^tp 4 , where we used that 
the dipole is of the order of l\ (electron cannot jump on a distance larger l\ due to 
exponential localization). The sum ^2 of the product of four random one-particle 
wave functions ip ~ 1/h runs over If sites for each electron so that ~ Vi an d 
U s ~ U/R 3 . The density of directly coupled two electron states is P2 ~ if/V and 
according to the Fermi golden rule the hopping rate is T 2 ~ U 2 p2 ~ U 2 lf/(VR 6 ). 
The mixing of two-electron levels takes place when n e = xl ~ ^iPi ~ ( r L^ 3 / r s) 2 > F 
where is the value of r s at the filling v = l/l 2 (one electron in a box of size l\ with 
Th = hU / (2^/nV)). For U ~ V this mixing takes place at R ~ l^ 3 ^> l\. In this sense 
the situation is qualitatively different from the case of short range interaction where 
the mixing is possible only for particles on a distance R < l\. For n e > 1 the pair 
of two electrons becomes delocalized in a way similar to the 3D Anderson transition. 
Indeed, the center of mass can move in the 2D plane and in addition electrons can 
rotate around a ring of radius R and width ij keeping their total energy E ~ U/R 
constant. Since R » l\ the dynamics is going in an effective 3D space where the 
Anderson transition to diffusion takes place at K e > 1. Formally, the diffusion of two 
electrons will be finally localized since the finite length of the ring makes the situation 
analogous to quasi-2D Anderson model (finite number of planes). However, the corre- 
sponding TIP localization length l c jumps from l c ~ l\ at n e < 1 to exponentially large 
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l c ~ /iexp(7rZ 1 K e ) 3> ii for K e > 1. In a similar way for M coupled 2D Anderson 
models the one-particle localization length l\ makes a jump near the 3D transition 
point (e.g. W/V = 16.5) from Zi ~ 1 to Zi ~ exp(g) where the conductance g ~ M. 
In the delocalized phase the diffusion rate or the conductivity in units e 2 /h can be 
estimated as D e ~ Zfr e ~ V^Kg/Zj. The above estimates show that the border between 
the insulating (localized at r s > r^ 3 ) and metallic (delocalized at r s < r^ 3 ) phases 
is: 

rs ~ ^ (2) 

The physical meaning of this result is very natural: larger r s corresponds to a lower 
density n s with larger distance between electrons and smaller two-body interaction 
between them. 

The first numerical studies of two polarized electrons in 2D Anderson model showed 
a transition in the level spacing statistics P(s) from the Poisson distribution Pp{s) 
near the ground state to the Wigner-Dyson statistics Pw( s ) at higher energies |27j 
in a way similar to the results of Shklovskii et al. |3Cj for the 3D Anderson model. 
In both cases the intermediate statistics P(s) are very close at the critical transition 
point (Fig. 2 in ^7j). The numerical studies were done in the range of parameters 
5 < W/V < 15, 0.1 < U/V < 2, 6 < L < 24, r s < 10; more details can be found in 
[ p7| p8| . Independent recent studies (3lJ of a similar model with Coulomb interaction 
at the center of the energy band also show a transition to P\y(s) as a function of 
W. Since even analytical formulas for the Coulomb matrix element between localized 
states are absent further investigations of this problem are still needed. 



3 Numerical studies of multi-electron problem in 2D 

For numerical studies of the Hamiltonian (1) we follow the approach developed in 
[ p7[ p9| and rewrite the hamiltonian matrix in the basis of one particle eigenstates 
(orbitals) at U = using the computed two-body matrix elements between the orbitals. 
Only a finite number M of low energy orbitals has been considered and the final many- 
body hamiltonian matrix was constructed on the basis of a pyramid rule for the one- 
particle orbital index mf. ^i=i m « — Si^i 1 * + M. This procedure allows to make 
an efficient reduction of the resulting matrix size N m comparing to mCn p without 
any serious modification of the properties of low energy states. To the maximum 
we used N m « 5000 corresponding to N p — 20, M — 42 (for N p = 30 we used up 
to M = 49). We checked that our results at low energy are not sensitive to the 
variation of M and N m (see the illustrations below). One of the main characteristics 
we extracted from the numerical studies is the level spacing statistics P{s) at a given 
total excitation energy E measured from the ground state. The disorder average 
has been performed over No = 5000 realization (for low E) and No = 1000 (for 
higher E). In this way, the total statistics used for P(s) varied from N s = 10 4 
at low E to Ns = 3 x 10 5 at high E. To study the variation of P(s) with the 
interaction or disorder strength it is convenient to use the parameter rj which is defined 
as rj = f°"(P(s)-P w (s))ds/ f*°(P P (s)-P w (s))ds, where s = 0.4729... is the smaller 
intersection point of Pp(s) and Pw(s). In this way 77 = 1 corresponds to Pp(s) and 
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Fig. 1 Dependence of 7/ on the rescaled total energy E/B for various numbers of particles 
N p = 3 (full triangle up), 6 (full diamond), 10 (•), 12 (triangle down), 20 (full triangle down) 
and 30 (*); W/V = 5, filling factor v « 1/32 and r s = 3.22, 10 < L < 31; B = W. 



T]=Q to Pw(s). We remind that for the one-particle 3D Anderson model the localized 
phase is characterized by r\ = 1, since localized (nonergodic) eigenstates do not feel 
each other, while in the metallic phase the eigenstates are extended and 77 = ppf . 

The results presented in [^9| show that at r s = 3.22 the parameter 77 evolves to 1 
for the total energy E < E c and for E > E c with the growth of the system size L 
at fixed filling v ~ 1/32. The critical energy E c and r\ c at which the transition takes 
place are E c w 0.25B, r] c w 0.56 (W/V = 10) and E c w 0.15S, 7y c w 0.33 (V^/V" = 7) 
with B = 4V^. For W/V = 15 no transition is found for E/B < 1 where 77 w 0.8 - 1.0. 
The new data for VF/F = 5 are presented in Fig. 1 where the transition is seen at 
E c w 0.1S, 0.19. At the critical point (E c ,r/ C ) the spectral statistics P(s) is 

independent of the number of particles and the system size reached in our simulations 
(see Fig. 3 in ]29]|). In Fig. 2 we show an illustration that at low energy the parameter 
77 is not sensitive to the variation of the total matrix size N m by more than 3 times 
that confirms the validity of our numerical approach. We note that not only the values 
of 77 are close for different N m but also the whole distributions P(s), as it is illustrated 
in Fig. 3 for r s = 5.76. 

The results of Fig. 1 are obtained at fixed parameter r s . The variation of the 
spectral statistics P{s) at low energy for the change of r s from 3.2 to 8.14 is shown 
in Fig. 3. These data demonstrate a clear approach to the Poisson distribution with 
growing r s . The same tendency is seen in the 77 dependence on energy shown at 
different r s in Fig. 4. These results are in agreement with the relation (2) and the 
argument given in the previous section according to which at large r s the distance 
between electrons grows and the two-body interaction between them drops. Since 
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Fig. 2 Dependence of rj on E/B for the case N p = 6 of Fig.l for different matrix sizes: 
N m = 1513, M = 26 (full diamond); N m = 498, M = 21 (square); here iV D = 5000. 




Fig. 3 Level statistics P{s) at VK/V = 7, L = 25 in the energy interval 0.1 < E/B < 0.14 for 
different r a : r a = 8.14, N p = 3, rj = 0.88 (•); r s = 5.76, N p = 6, r? = 0.77, M = 27, iV m = 1844 
(x) (compare with (+) at smaller basis M = 24, N m = 996, 77 = 0.78); r s = 3.15, N p = 20, 
T) = 0.52 (o); total statistics Ns > 2.5 x 10 4 . Full lines show the Poisson distribution and the 
Wigner surmise. 
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Fig. 4 Dependence of r\ on E/B for different r a at L = 25: r a — 8.14, N p = 3 (full triangle 
up); r s = 5.76, iV p = 6 (full diamond); r s = 4.07, N p = 12 (triangle down); r s = 3.15, 
N p — 20 (full triangle down). 

only the interaction between two particles is able to mix many-body levels, the 
electron dynamics becomes nonergodic at large r s that leads to the Poisson distribu- 
tion. 

At the same time even at the optimal values of r s w 3 the ergodicity and level 
mixing appear only at total energy E > E c while for E < E c including the ground 
state the ergodicity is absent in the range of parameters we studied. Formally, one can 
object against the statement that the value rj = 1 implies the localization of electrons 
in space giving as a counterexample the case of noninteracting delocalized electrons 
in the metallic phase where also rj — 1 for the multi-electron spectrum. However, 
in the system we study all one-particle states are exponentially localized and it is 
very difficult to imagine that the interaction gives derealization of charge, makes 
the system metallic and does not introduce complete ergodicity and level mixing, 
which should give 77 = instead of 77 = 1 found for E < E c . Due to this reason 
we consider that the appearance of Poisson statistics near the ground state is the 
direct evidence of localization at E < E c . At the same time the emergence of random 
matrix statistics at E > E c implies ergodicity and derealization for many electron 
states. This derealization apparently appears in a way similar to the two electron 
derealization discussed in the previous section but in addition the interaction between 
larger number of particles decreases the derealization border comparing to the TIP 
case discussed in [ p7| p8| . An interesting point is that the transition at E — E c takes 
place at zero temperature. Indeed, according to the data of Fig. 1 at the transition the 
energy per particle e c = E c /N p «5x 10~ 3 B is much smaller than the Fermi energy 
ep — 0.1B showing that the data are close to the thermodynamic limit with T = 0. 
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4 Conclusion 

Our results show that the Coulomb interaction between two electrons in excited states 
leads to their derealization for 1 < r s < r^ 3 while for r s > r^ 3 they remain localized. 
The transition between these phases is similar to the Anderson transition in an effective 
dimension d e f / = 3. The numerical studies of the spectral statistics for many polarized 
electrons in the 2D Anderson model show that for 3 < r s < 9 and 5 < W/V < 
15 the ground state is nonergodic (localized) and is characterized by the Poisson 
statistics for the total energy E < E c . However, the transition to quantum ergodicity 
and the Wigner-Dyson statistics takes place at a fixed total energy E c independent 
of system size (for r s w 3.2, 5 < W/V < 10 and fixed filling v w 1/32). This 
implies a derealization at zero temperature T. At the critical point E c the parameter 
rje increases with the disorder strength W: r\ c = 0.19 (W/V = 5); 0.33 (W/V — 
7); 0.56 (W/V = 10) and the critical statistics approaches to the Poisson limit. In 
a certain sence the situation is similar to the Anderson transition in high dimensions 
d > 3 discussed in Q where also r\ c is growing with d. In analogy with this result and 
the case of two electrons in 2D we make a conjecture that in the Hamiltonian (1) the 
transition at E c is similar to a transition in some effective dimension 3 < d e ff < 2N p . 
This d e ff is growing with the disorder strength W. 

The interaction induced ergodicity at T = is in favor of the metal-insulator 
transition observed experimentally, especially in the view of recently observed metallic 
behavior for polarized electrons [fl5| . However, our data are not sufficient to determine 
the behavior of resistivity on temperature in the ergodic phase at E > E c . Therefore, it 
is not excluded that at strong disorder this ergodic phase will show a resistivity growth 
at low T. In this case one can suppose that the ergodicity induced by the Coulomb 
interaction is responsible for the phononless VRH conductivity as it was argued in 
[ p"9[ [T^] and indicated by the experiments [Tt]]. More detailed investigations are 
required to understand the properties of the Coulomb ergodic phase at E > E c . 

We acknowledge the IDRIS at Orsay for the allocation of the CPU time on supercomputers. 
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